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Abstract 

We develop a semiclassical theory for the dynamics of electrons in a magnetic 
Bloch band, where the Berry phase plays an important role. This theory, 
together with the Boltzmann equation, provides a framework for studying 
transport problems in high magnetic fields. We also derive an Onsager-like 
formula for the quantization of cyclotron orbits, and we find a connection 
between the number of orbits and Hall conductivity. This connection is em- 
ployed to explain the clustering structure of the Hofstadter spectrum. The 
advantage of this theory is its generality and conceptual simplicity. 

PACS numbers: 72.10.-d 72.15.Gd 73.20.Dx 
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The theory of semiclassical dynamics of Bloch electrons in a weak electromagnetic field 
plays a fundamental role in our understanding of electronic spectral and transport properties 
in metals and semiconductors. The basic ingredients of this theory are the following pair of 
equations: 

r=^, ^k=-eE-efxB, (1) 

where £ n (k) is the energy for the n-th band, and E and B are the electric and magnetic 
fields II . The validity of these relations depends on the absence of interband tunneling, 
therefore (1) holds if the external fields are sufficiently weak. 

A natural question is how should these semiclassical equations be modified for a magnetic 
Bloch band (MBB). Such a band is obtained when an electron is subject simultaneously to 
a periodic potential and a magnetic field (not necessarily weak), such that the magnetic flux 
per unit cell of the periodic potential (plaquette) is a rational multiple of the flux quantum 
h/e. For example, in a tight binding model Hofstadter showed that a Bloch band is broken 
into q subbands if the rational number is p/q ||. At the opposite limit, when the magnetic 
field is much stronger than the periodic potential, a Landau level is broadened and split into 
p subbands ||. These subbands at both limits are manifestations of the magnetic Bloch 
bands. 

It is difficult to resolve these MBBs in a naturally occuring solid. For example, with 
a lattice constant a = 5 A and a magnetic field B = 3 Tesla, the value of p/q is of order 
10 -2 . However, p/q can be a significant fraction of unity if we are using an artificial lattice 
with a much larger period. This is realized, for example, by imposing an optical interference 
pattern on top of a heterostructure to create a grid potential on the two dimensional electron 
gas beneath. Recent experiments have demonstrated that disorder can be reduced to such a 
degree that the effect of magnetic bands emerges in the transport properties of the electron 
gas [|. 

Our goal is to derive the counterpart equations to Eq. (1) for a MBB, and use them to 
explore the dynamics of electrons under a weak electromagnetic perturbation. To simplify 
the discussion, assume the electrons are confined in a two-dimensional periodic potential, 
with a magnetic field perpendicular to it. Because the vector potential for a constant 
magnetic field is not periodic, the Hamiltonian H does not commute with usual translation 
operators. However, we can define magnetic translation operators T(R) that commute with 
H H . Analogous to the Bloch states, we require the eigenstates of H to satisfy the relation 
T(R)$ n (k) = exp(zk-R)\I/„(k) ||. This is not well-defined in general, because the magnetic 
translation operators for different displacements, Ri and R2, do not commute unless there 
is an integer number of flux quanta in the area |Ri x R 2 |. Therefore, if = p/q, we have 
to choose a unit cell consisting of q plaquettes. Correspondingly, the area of the magnetic 
Brillouin zone is reduced by a factor of q. Furthermore, because of the magnetic translation 
symmetry, the energy spectrum is exactly g-fold degenerate. 

Electric perturbation and transport - - If we write a magnetic Bloch state in the form 
\I / n (ko) = exp(ik ■ r)u n (k ), then u n is modified by a weak and homogeneous electric field 
into 
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where k = k(i) = ko — eEt/h, and we have used a time dependent vector potential for the 
electric field. The above result is obtained by adiabatic perturbation theory, which is valid 
for weak fields. The average velocity in such a state can be easily evaluated as: 

r = ^^--kxzQ n (k), hk=-eE, (3) 

where z is the unit vector along the direction of the magnetic field. The second term in the 
expression for r comes from the first order non-adiabatic correction in the wave function, 
with § 

rwi \ . ( ,du n du n du n du n \ 

The equations in Eq. (3) are the new set of semiclassical equations. Notice that the usual 
Lorentz force term for k is absent because the magnetic field has already been included in the 
band structure. On the other hand, the velocity has an extra term involving fl n (k), which 
will be called the "curvature" of the Berry phase, because its integral over an area bounded 
by a path C in k-space is the Berry phase T n (C) ||. In physical terms, fi n (k) describes 
the contribution of state \l/ n (k) to the Hall conductivity in the absence of scattering. The 
derivation of Eq. (3) is based on a homogeneous field; nevertheless, it should still be valid 
when the field is slowly varying in space and time. The generalization of Eq. (3) to higher 
dimensions is straightfoward. In that case there will be more than one component of the 
curvature. 

The combination of Eq. (3) with the Boltzmann equation, 



dr <9k 



coll 



offers a general framework for semiclassical transport in a MBB. The right hand side is the 
collision term due to impurity scatterings, etc. It has to be cautioned that the Boltzmann 
equation is valid only when the scattering broadening of a MBB is small compared with 
its bandwidth. This does not pose essential difficulty when p/q is a simple fraction. In the 
more general situation of a large or infinite q (irrational <ft), we have to modify our approach 
in the following way: We divide the total magnetic field B into -Bo and SB, where B relates 
to the band structure not destroyed by disorder, and SB is a small perturbation. Then the 
semiclassical dynamics in the MBBs of -Bo, driven by E and SB (see Eqs. (7) and (9) in the 
next section), will be employed in the Boltzmann equation. Under such a circumstance, the 
scattering broadening is only required to be smaller than the bandwidth for B . 

To demonstrate the use of Eq. (3) in transport problems involving an electric perturba- 
tion, we consider a homogeneous system in which / depends on k only, and use the relaxation 
time approximation. The current to first order in E is: 
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E x z*j- J ^f Q n (k) (6) 




where r(k) is the relaxation time and /o is the unperturbed Fermi-Dirac distribution. The 
first term is new and is due to the Berry phase curvature, which is nonzero in general. In 
fact, in the simple case of a filled band (fo = 1), for which the second term is zero, this term 
reduces to the topological Chern number discovered by Thouless et al J7J. The second term 
is the usual Boltzmann transport formula . It was used for the calculation of longitudinal 
conductivity of magnetic bands ||. Our theory justifies this usage because the Berry phase 
term only contributes to the Hall conductivity. 

Compared to the brute force, all purpose Kubo formula approach, the semiclassical 
dynamics, in conjunction with the Boltzmann transport theory, offers a simple and intuitive 
picture of the behavior of the physical system. The semiclassical dynamics also provides a 
useful tool in problems with spatially varying and/or time dependent fields, where quantum 
mechanical calculations are usually very involved. A detailed study will appear in a separate 
publication |IU[ . 

Magnetic perturbation and hyperorbits — The remaining part of this letter will focus on 
how the MBBs for a given magnetic field B are perturbed by adding SB. In this case, the 
equations for the semiclassical dynamics become 

BE (k) 

r= ^ ; -kxfO B (k), hk = -er x SBz. (7) 
This result can be derived, for instance, by considering a wave packet in a MBB and studying 



how its center of mass moves in r-space and k-space ||10|| . Notice that the wave vector k, 
which is a good quantum number for Bo, is no longer conserved in the presence of SB, even 
when both B and SB are uniform. 

After r is eliminated by combining both equations in Eq. (7), the equation for k takes 
the following form: 

• _ dS n /dk x zSBe/h 

1 + Q n {k)SBe/h ' 1 j 

It is not difficult to see that k moves along a constant energy contour in the magnetic 
band structure. The presence of fi n (k) changes the speed of motion, but it does not alter 
the shape of the orbit for a given energy. The cyclotron orbit in r-space can be derived 
from f = kx z(h/eSB), which shows that the r-orbit is simply the k-orbit rotated by n/2 
and scaled by the factor Ti/eSB. Such "hyperorbits" were introduced by Pippard [[nj. We 
emphasize that, the existence of hyperorbits is a quantum effect, and cannot be explained 
classically. One possible way to detect them is by using an electron focusing device with a 
configuration similar to a mass spectrometer [JlIJ . In order to have a successful observation, 



the hyperorbit has to be within the ballistic range of electron transport. 

Analogous to ordinary cyclotron orbits, these hyperorbits will drift in an external electric 
field. By adding a — eE term to the second equation in Eq. (7), we obtain 
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h ■ E x z , . 

The time average of the first term for a closed orbit is zero, while the second term describes 
the drifting of the hyperorbit that results in a Hall current. This will be used later to 
calculate the Hall conductivity for magnetic subbands in the Hofstadter spectrum. 

Quantization of hyperorbits — The dynamical equation for k in Eq. (8) can be cast into the 
Lagrangian formulation 

h 2 

L(k, k) = ^(hk 2 - hh) - £ n (k) + HA n ■ k, (10) 

where A n is the "vector potential" for the Berry phase that satisfies V x A n (k) = Q n (k)z. 
Apart from an unimportant constant, the propagator for a completed closed orbit C in 
period T is given by exp(i/h J T Ldt). For a semiclassical orbit, the amplitudes for paths 
that circle different times must add constructively. This leads to the following quantization 
rule for the area of a hyperorbit in the n-th MBB , 

-I (kxdk). Z = 2n U+- _ _ , (11) 



where m is a non-negative integer, and T n (C m ) is the Berry phase for orbit C m [14 . 

By using the constraint that the area of the outer-most orbit be smaller than the area of 
the first magnetic Brillouin zone of size (2ir/a) 2 /q, we found the number of quantized orbits 
to be the integer part of l/(q5(f>) + T n {C max ) /2n + 1/2, where 5(j) = 5Ba 2 e/h. T n (C max )/27i 
can be replaced by the Hall conductivity a n (in units of e 2 /h), because the Berry phase for 
the orbit C max is very close to 2na n . For an integer value of l/(g<50), we then have 

1 



number of orbits 



+ O r 



(12) 



These orbits will be broadened into subbands by tunneling to orbits with the same energy 
in other magnetic Brillouin zones. Therefore, this naive-looking formula relates the Hall 
conductivity a n of a parent band to the number of daughter subbands under a perturbation 
5<p. It is crucial in understanding the clustering pattern for the Hofstadter spectrum. 

The Hofstadter spectrum — Consider a Bloch band subject to a magnetic flux that can be 
expanded as 

• (13) 



fi + 



Its r-th order approximation will be written as r = p r /q r - At the first order, the Bloch 
band is broken into f\ subbands; each subband is further fragmented by an extra magnetic 
field at the second order, and so on. As will be shown below, the new semiclassical dynamics 
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offers a clear and intuitive picture about how each subband (parent) in the r-th order should 
split into daughter subbands fli7| . 

Our major findings are summarized below: (a) Firstly, it is important to distinguish 
between "closed" and "open" subbands. We define a "closed" subband to be a subband 
broadened from a closed hyperorbit; similarly an "open" subband is derived from an open 
orbit. For a square or a triangular lattice, all subbands except one for every parent band 
are closed, (b) The closed subbands at the same order all have the same Hall conductivity 

a r = (-ir-Vi, (14) 

where the subscript in a refers to the order, not to the band index of the subband. If a 
parent band has only one open daughter band (eg., for a square or a triangular lattice), 
the Hall conductivity for this open subband is (— l) r_1 g r _i + (— l) r q r - (c) Because of the 
difference in the Hall conductivities, a closed band will break into f r+ % subbands at the next 
order, while an open band will break into f r+ \ + 1 subbands. 

We briefly describe the derivation for the Hall conductivity of a closed subband. It is 
determined by the response of the corresponding hyperorbit under an electric field. Using 
Eq. (9), we know that the drifting velocity of a closed orbit is (r) = — E x z/8B. It follows 
that the Hall conductivity for a closed subband at the r-th order is o r = ep r /5B r _i, where 
p r is the electron density per unit area and 5B r _i = M0. r _i/(ea 2 ). Since the electrons are 
equally distributed among the q r subbands at the same order, p r is equal to l/q r times 
the electron density of the original Bloch band. Eq. (14) is obtained after the identity 
p r /q r — Pr-i/q r -i = (— l) r_1 /(9r9r-i) is used to evaluate 5(p r -i. The Hall conductivity for 
an open subband can be figured out quite easily by using the sum rule: cr par e n t = ^daughter 
H 0- 



For a square lattice, an open daughter band is always located at the center of a parent 
band; therefore, we know which subband the conductivity cr open belongs to. The Hall con- 
ductivity distribution obtained this way is exactly the same as that which is obtained by 
using the Diophantine equation with some subsidiary constraints |7j]. We should emphasize 
that this is the first time the (seemingly) erratic behavior of the Hall conductivities for the 
Hofstadter spectrum is given a clear and direct physical meaning. 

With the help of Eq. (12), we can determine how a parent band is splitted by 5(f). Sub- 
stituting 5(p r = (— l) r / (q r+ xq r ) into Eq. (12), and taking into account the g r -fold degeneracy 
for the r-th order magnetic Brillouin zone, we then have the number of daughter bands for 
an r-th order parent band 

V r = \(-l) r q r+1 + a r \/q r . (15) 

In conjunction with Eq. (14), we obtain T> r = f r+ i for a closed band. Similarly we have 
T> r+ i = f r+1 + 1 for an open band. Azbel conjectured that there are f r+ i daughter subbands 



for every parent band [ 17j • It is clear from our calculation that his conjecture is correct only 
for a closed parent band. An expression similar to Eq. (15) has been obtained by Wilkinson, 
but his evaluation of such an expression required external input for the Hall conductivity 

Various values of have been used to check the predicted splitting of subbands with the 
actual Hofstadter spectrum, and the results are found to agree very well. One exception 
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occurs when there is a degeneracy among the subbands. In this case the individual oy, 
as well as T> ri cannot be determined uniquely. For example, when = 1/4, the central 
two subbands are degenerate for a square lattice and the Hall conductivities can be either 
(1, 1, —3, 1) or (1, —3, 1, 1) (q = 1). One of these two subbands will have this —3 (cr open ) 
when the degeneracy is lifted by next nearest neighbor coupling |L9| . As expected, we found 
the subband corresponding to a hyperorbit that is closer to the nesting energy contour 
acquires this <7 open . 

For a triangular lattice, the open orbit is near the zone boundary (exact location requires 
some calculation). Therefore, the distribution of Hall conductivities and band splitting are 
no longer symmetric in energy. This explains the asymmetry in the spectrum generated by 
numerical calculation [19|. For a lattice without 3-fold or 4-fold symmetry, more than one 
open orbit may exist in a range of energy. In this case, the total Hall conductivities for open 
orbits can be figured out by the sum rule. However, further analysis is required to get the 
detailed distribution within them. 

In summary, we have demonstrated the new semiclassical dynamics for magnetic Bloch 
bands and its application to a variety of phenomena involving strong magnetic fields. It can 
be used to calculate the transport properties, to obtain the quantization rule for hyperorbits, 
and to get the Hall conductivity for a magnetic subband. We also showed that the complex 
structure of the Hofstadter spectrum can be explained as a logically consistent part of this 
theory. 
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